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fortunate situations, overlapping patterns of two or more incomplete series may give 
complete sequence information. Ideally, one would prefer to observe complementary 
information from series of N-terminal and C-terminal fragment ions to bolster 
confidence in the analysis. 

5 [0059] Consider the MALDI-PSD mass spectrum shown in Figure 8 as an 

unknown. It can be assumed at the outset that the major peak at m/z 574 Da 
represents the protonated molecule, which was the precursor ion selected for PSD. 
The protonated molecules fragment during the PSD process and degrade into 
fragment ions represented by the peaks at lower m/z shown in Figure 8. The 

10 procedure for analysis or data interpretation consists of merely examining the mass 
difference between each of the fragment ion peaks and the peak representing the 
protonated molecule. The goal is to find a fragment ion peak that differs in mass 
from the protonated molecule peak by either a residue mass or a residue mass plus 
water. A peak at a mass-to-charge (m/z) value that differs from the protonated 

15 molecule peak by the mass of a residue mass plus water corresponds to the amino 
acid that was located at the C-terminus of the original peptide. 

[0060] In considering the mass spectrum in Figure 8, the peak at m/z 425 
corresponds to the loss of 149 Da, which corresponds to the mass residue of 
methionine and water; i.e., 149 = 131 (the residue mass in methionine) + 18. The 

20 observation of the peak at m/z 425 suggest that methionine was expelled from the C- 
terminus to form a b ion at m/z 425. Having recognized the largest b ion in the mass 
spectrum by its peak at m/z 425, the goal is to identify the next smaller b ion, namely 
one that also has lost the second amino acid residue from the C-terminus to form the 
second b ion. The peak at m/z 278 represents a mass difference from 425 that 

25 corresponds to only the residue mass of an amino acid. That is, the mass difference 
between 425 and 278 is 147 Da, which is the residue mass of phenylalanine. This 
suggests that the second amino acid in the sequence from the amino terminus is 
phenylalanine: -Phe-Met. 

[0061 ] Continuing with the b series, the peak at m/z 221 is 57u lower than 
30 278. This mass difference corresponds to the residue mass of glycine: -Gly-Phe- 
Met. This observation suggests that the third amino acid from the C-terminus is a 
glycine. residue. Finally, it can be seen that there are no other peaks observed below 
m/z 221 that correspond to the loss of a residue mass from 221; the search for 
further b ions come to a halt. At this point, the peak represents the intact molecule 
35 (the peak at m/z 574 which is the protonated peptide) and an examination of 
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fragment ion peaks that differing in mass from 574 by exactly a residue mass to try to 
recognize the beginning of the y series of ions, is pursued. Upon reconsideration of 
this mass spectrum, one notices that the peak at m/z 41 1 is 1 63 lower than the peak 
, at m/z 574. This corresponds to the residue mass of tyrosine. Residue loss from the 
5 protonated molecule suggests that tyrosine is the amino acid located at the amino 
terminus of the peptide. 

[0062] With the focus now on m/z 41 1 as a reference point, notice that 
the peak at m/z 354 is 57 Da lower, a mass which corresponds to the residue of 
glycine. The observation of peak at m/z 354 suggests that the second amino acid 

10 from the N-terminus is a glycine residue, giving the partial sequence: Tyr-Gly-. 
Using the peak at m/z 354 as a reference point, the peak at m/z 297 is 57 Da lower 
than m/z 354 and suggests the expulsion of another residue of glycine in the 
formation of this ion. These observations suggest that the third amino acid residue is 
also a glycine residue in sequence from the N-terminus. Continuing with the peak at 

15 m/z 297 as a reference point, the peak at m/z 150 is 147 Da lower, suggesting 
expulsion of a phenylalanine residue or that phenylalanine is the fourth amino acid 
from the N-terminus, giving the partial sequence: Try-Gly-Gly-Phe-. Since there are 
no other peaks at lower m/z that differ from m/z 1 50 by a residue mass of an amino 
acid, recognition of the sequence comes to a halt. 

20 [0063] Having suggested possibilities for amino acid sequences from 

either terminus of the peptide, these two suggested partial sequences are overlaid to 
compose a complete sequence. For example, the data reviewed above suggest that 
the sequence at the C-terminus is: -Gly-Phe-Met; whereas another series of data 
suggested that the series of amino acids at amino terminus was: Try-Gly-Gly-Phe-. 

25 Note that each of these two partial sequences shows the C-terminus of a glycine 
residue connected to the amino terminus of a phenylalanine residue. These two 
residues must be a redundant observation in the two sequences, and they can be 
overlaid at that point. This would give an overall sequence starting from the amino 
terminus of. Try-Gly-Gly-Phe-Met. . 

30 [0064] If the residue masses of the postulated amino acids in the 

complete sequence are summed, 555 Da will be the obtained value. Add 18 Da to 
this sum for the hydrogen at the amino terminus and the hydroxyl at the C-terminus. 
Then add 1 Da for the proton on the protonated molecule. This will give a total of 
574 Da for the expected mass of the protonatetd molecule, in perfect agreement with 
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the observed peak at m/z 574, giving credence to the suggested amino acid 
composition implicit in the sequence Y-G-G-F-M. 

[0065] The protein sequencing can be conducted by use of the present 
invention with or without use of a matrix. The use of ah ultrafast, femtosecond laser 
5 is envisioned to minimize any destructive burning of the specimen, thereby potentially 
rendering use of an expensive and time-consuming matrix as unnecessary. Without 
a matrix (herein, also known as having "isolated molecules"), the identification and 
sequencing of the protein is simplified since the matrix characteristics do not have to 
be accounted for and filtered out of the calculations. A femtosecond laser in the 

10 range of approximately 20 femtosecond duration pulses allows for localization of the 
energy based on the speed of the pulse and the ability to quickly shape the phase 
and amplitude modulation of the pulse. Furthermore, the specimen fragmentation is 
primarily due to laser cleavage rather than enzyme or chemical cleaving. This is 
ideally suitable for insoluble proteins, glycocylated proteins which have been linked to 

15 cancer, (including the selective cleavage of the associated oligosacharides) direct 
protein analysis from silicate substrates, direct analysis of' PAGE gels, direct 
sampling of membrane proteins from intact cells and bacteria, the direct sampling of 
. genetically modified agricultural produce (such as grains), and even human matter 
such as hair, fingernails and fingerprints. 

20 [0066] The personal computer employs a method and software for protein 

sequencing as follows. The foundation of this method is based on the fact that there 
are only 20 amino acids and that their masses are well known. First, the computer 
determines the molecular weight of the intact proteins specimen. This requires the 
generation of a single high-mass peak and minimization of the low weight 

25 background. Secondly, the computer automatically finds peaks that are an integer 
number of amino acids smaller than the parent protein; a laser beam pulse shape 
that causes some fragmentation can be employed. Thirdly, this procedure is 
continued from high to low masses. Finally, confirmation of results can be 
automatically obtained by a statistical optimization method (such as that previously 

30 described for the MALDI process) that attempts to optimize a given mass; the 
success of this optimization will depend on whether that fragment of the protein has 
an integer number of amino acids. Automatic adjustment for the N or C terminus is 
also automatically adjusted for by the computer as previously explained. Alternately, 
each single amino acid could be separately searched for. Thus, the present 

35 invention control system and apparatus is ideally suited for analyzing, identifying, 
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sequencing and severing complex multi-molecular specimens in a highly automated 
manner. 

Selective Bond Cleavage 

[0067] The ultra-fast laser of the present invention is used to enhance in- 
5 source photochemistry and fragmentation, however, random fragmentation would not 
be as useful as selective bond cleavage. Furthermore, selective peptide bond 
cleavage would be ideal for protein sequencing. Cleavage of amino acid side chains 
may be of value for de novo sequencing because it would allow a determination of 
the presence or absence of certain amino acids. Similarly, selective cleavage of 

10 phosphate groups, oligosacharides and other post-translational modifications would 
be equally valuable. The ideal, of course, would be to achieve peptide bond 
cleavage without loss of side chains or other appended 'groups. This would allow, for 
example, to determine phosphorylation sites. 

[0068] It is envisioned that selective bond cleavage can be realized when 

15 using shaped pulses that are capable of localizing the energy in a time scale that is 
short enough to prevent total energy randomization. For example, the protonated 
molecule of bradykinin potentiator C, as produced by MALDI, fragments poorly during 
PSD, and does not produce a suitable spectrum from which one could deduce the 
amino acid sequence. Thus, this 11 -residue peptide is ideal for this application. 

20 Selective laser bond cleavage may have additional application as a synthetic route to 
specific products. 

Fixed Shape Pulse Shaping Apparatus 

[0069] A second preferred embodiment of a laser excitation apparatus 
121 is showing in Figure 12. Apparatus 121 includes a femtosecond laser 123, an 

25 upstream prism -125, a downstream prism 133, a pulse shaping mirror 129 at the 
Fourier plane, an offset mirror 131, and a receiver or targeted specimen 135. 
Upstream prism 125 initially acts to disperse the colors of the emitted laser beam 
pulse while downstream prism 133 serves to redirect this dispersed laser beam pulse 
toward pulse shaping mirror 129. Pulse shaping mirror 129 has a predetermined or 

30 fixed pulse shaping surface, such as a set of 800 nm wavelength sine curves, formed 
or machined therein. The physical characteristic or shape of the actual pulse 
shaping surface is predetermined through optimization experimentation for the 
intended use and intended laser beam input by use of a learning program such as 
that disclosed in the first preferred embodiment. After the desired mirror surface 

35 shape is known for the intended use, the less expensive, fixed shape mirror 129, or 



WO 02/061799 



PCT/US02/02548 



an alternate fixed pulse shaping optic, can be employed to reduce equipment costs 
for actual production systems. Also, the computer and optimization program are not 
required for these types of known set up and known specimen applications after the 
initial determination is conducted. 

5 [0070] The passive pulse shaping mirror 129 thereby reshapes the laser 

beam pulse shape characteristic, reflects it back through the same prisms in reverse 
order, and in an offset or time-delayed manner. The position and orientation of mirror 
129 alters the final characteristics of the shaped pulse. Thus, a computer controlled 
automatic actuator can move or change the position or orientation of mirror 129, 
10 based on an algorithm, such as shaping the pulse with the cosine portion of the sine 
wave shape. Offset mirror 131 subsequently reflects the shaped laser beam toward 
the receiver, which can be a mass spectrometer, fiber optic sensor/switch, or a 
targeted specimen. 

[0071 ] It is alternately envisioned that an in-line optical apparatus can be 
1 5 used, such as that disclosed with the first preferred embodiment, however, the pulse 
shaper at the Fourier plane would be replaced by a passive mask having a 
transmissive optic with a predetermined coefficient of refraction, or a polarizing-type 
sine mask on a transparent substrate. Also, a polymer-doped glass or blend of 
polymer sheets that are capable of retarding the phase of the laser beam pulse wave 
20 or otherwise varying a wavelength, timing or shaping characteristic of same can be 
employed. 

[0072] Alternately, certain optics can be used such as a backside coated, 
chirped mirror having multiple dichroic layers, which wpuld be satisfactory for pulse 
shaping without dispersive" optics and without the need for a Fourier plane. An 

25 acceptable chirped mirror is disclosed in Matuschek, et al, "Back-side-coated Chirped 
Mirrors with Ultra-smooth Broadband Dispersion Characteristics," Applied Physics B, 
pp. 509-522 (2000). A negative dispersion mirror from CVI Laser Corp., part no. 
TNM2-735-835-1037 is another suitable example. A rotatable wheel having multiple 
different chirped mirrors, each with specific pulse shaping characteristics, can also be 

30 used to provide a discrete number of predetermined shaped pulses. 
Optical Coherence Tomography 

[0073] A third preferred embodiment of the present invention uses an 
apparatus 221 for laser excitation or ionization with Optical Coherence Tomography 
("OCT). In general, Figure 17 illustrates the OCT application of apparatus 221 
35 wherein there is a femtosecond laser 223, a laser beam shaper 229, a human or 
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animal tissue specimen 241 , an optical gate 251 and an image 253. Laser 223 emits 
a laser beam pulse shorter than 1 picosecond. Shaper 229 is made of three parts; 
two dispersive elements 255 which sandwich a phase mask element 257. Shaper 
229 essentially prevents multiphoton excitation which can damage the person's or 
5 animaFs DNA, as will be discussed in more detail as follows. An unshaped laser 
beam pulse is used to gate the ballistic photons to render the image for tomography 
use. Optical gating can be accomplished by up-conversion in a frequency doubling 
crystal or with a kerr-gate in liquid carbon disulphide. The construction of apparatus 
221 as illustrated supposes transmission imaging; the same end result can 

10 alternately be accomplished with back scattered imaging. Image 253 could be 
viewed like an x-ray-type image of the internal organs of the human or animal 
specimen but without harmful three photon exposure. The use of the shaped pulse 
in OCT provides for an increase in laser intensity for better imaging while preventing 
the damaging effects caused by multiphoton excitation of healthy tissue. 

15 Photodvnamic Therapy 

[0074] A fourth preferred embodiment of the present invention uses an 
apparatus also shown as 221 for laser excitation or ionization with photodynamic 
therapy ("PDT"). . In general, Figure 17 also illustrates the PDT application of 
apparatus 221, but optical gate 251 and image 253'are not being required. Shaper 

20 229 allows two photon excitation but essentially prevents three photon excitation. 
Shaper 229 enhances the laser induced activity of a therapeutic agent which 
prevents damage of healthy tissue. Use of laser beam pulse shaping of the present 
invention should provide superior control and results for PDT applications as 
compared to those practically possible with conventional methods as disclosed, for 

25 example, in U.S. Patent No. 6,042,603 entitled "Method for Improved Selectivity in 
Photo-Activation of Molecular Agents" which issued to Fisher et al. on March 28, 
2000, and is incorporated by reference herein. Alternately, the pulse shaper can be 
tuned to target cancerous cells for multiphoton gene therapy or destruction, with or 
without the presence of a therapeutic agent, without damaging healthy tissue. 

30 Contro l .tf Nonlinear Optical Processes 

[0075] As applied to all of the applications herein, selective control of one 
and multiphoton processes in large molecules, including proteins, is possible using 
simple pulse shaping. The results show an extraordinary level of control that is robust 
and sample independent, with contrast ratios near two orders of magnitude (clearly 

35 visible with the naked eye). Such large contrast ratios allow for more precise 
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cancellation control of undesired photons and other laser beam characteristics, such 
that nonlinear transitions induced by each pulse are controlled. Because simple 
phase functions can be incorporated into a passive optical component such as mirror 
129 (see Figure 12), these applications do not require the complexity and expense of 
5 computer controlled pulse shapers after initial set up, although systems such as in 
Figure 1 can still be employed. 

[0076] The underlying concept of the apparatus and associated method 
are shown in Figures 13A - 13C. Multiphoton transitions are optimized when the 
central bandwidth of the laser pulse u) 0 , is some fraction (half for two-photons, a third 

10 for three-photons, etc.) of the total energy of the transition as illustrated in Figures 
13A and 13C. For ultrafast pulses, when the bandwidth is large, different frequency 
components (<d 0 ± Q) of the pulse can interfere, thereby reducing the probability for 
multiphoton excitation. Referring to Figure 13B, the spectrum of the ultrafast laser 
pulse with amplitude A(Q.) is plotted as a function of detuning from the central 

1 5 frequency. A phase mask <p(Q) can be imprinted on the pulse such that the phase of 
each frequency component Q. acquires a specific value. The effect of pulse shaping 
on the probability for two-photon absorption ("2 PA") can be calculated as follows: 



Pip*?* 



[4] 



DO 

Ja(&) A(-&) exp[i{<p(Q) + <p(-Q) } ]dQ 



20 and for three-photon absorption ("3PA"), a similar formula can be derived as follows: 
[5] 



J jA(Q 1 )A(Q 2 )A(-£2 l -Q 3 )exp[i{<p(Q 1 ) + q)(a a ) + q)(-G 1 -Q 2 )}]d£2,d£2 : 



where amplitudes and phases are introduced for two different detuning values Q and 
Q 2 , as shown in Figure 13C. One photon transitions are not affected by the phase of 

25 the pulses, however, one photon processes are difficult to selectively achieve at high 
photon flux due to the onset of multiphoton processes. 

[0077] A schematic representation of two photon and three photon 
induced fluorescence is illustrated in Figures 13A and 13B, respectively. The vertical 
arrows represent ultrafast pulses that induce the two and three photon transitions. 

30 Because of their broad bandwidth, ultrafast pulses contain photons that are detuned 
from the central wavelength © 0 by an amount Q. Referring again to Figure 13C, 

Ultrafast laser Dulses are sVia.n«rt nsinn a sine fiinr!ttr>n nhaeo maoU 9m eo *h* m.ie* 
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spectrum underlying the dashed curve while the structures of the chromophores are 
. . also shown. 

Example 1 

[0078] The experiments in all of the following, example were carried out 
5 using an amplified titanium sapphire laser producing 50 fs pulses. The pulses were 
shaped using a spatial light modulator ("SLM") at the Fourier plane of a zero- 
dispersion two grating arrangement. The two independent modulator plates, based 
on liquid crystal technology in the SLM (128 pixels each), were calibrated so that only 
phase delays were introduced without changes to the output spectrum, intensity, and 

10 polarization. The shaped pulses centered at 809 nm were characterized by second 
harmonic generation frequency resolved optical gating. When all phases were set to 
zero, laser pulses were near transform limited. Unless indicated otherwise, 
measurements were made with pulse energies of 0.4 jd/pulse at the sample. 
Experiments were carried out by setting the phase function equal to a sinusoid, as 

15 shown in Figure 13B, in the 779-839 nm spectral range. Emission from one photon 
or multiphoton induced processes from various samples was measured as a function 
of 5, the phase shift of the mask across the spectrum. The maximum phase 
advancement or retardation was 1 .5jc. 

[0079] Equations 4 and 5 can be used to calculate the expected signal for 

20 two and three photon processes as a function of 8. These calculations are graphed 
in Figures 14A - 14H for sinusoidal phase functions having a half (Figures 14A and 
14B) or a full (Figures 14C and 14D) period across the entire phase mask. The 
calculated probability for two photon and three photon transitions peaks at half 
integer values of % in Figures 14A and 14B, while the calculated probability for two 

25 photon and three photon transitions peaks at integer values of n in Figures 14C and 
14D. The shape of the phase function, where maxima and minima in the probability 
are achieved, is indicated as inserts. 

[0080] Experimental data were obtained with the phase functions used for 
the calculations in Figures 14A - 14D. In these experiments, the two and three 

30 photon emission from large organic molecules is detected as a function of 5. 
Although the model described by equations 4 and 5 assumes two level systems, 
Figures 14E -14H experimentally demonstrate that this principle can be applied to 
complex systems having a manifold of vibrational states broadened by the presence 
of a solvent. It is noteworthy that the peaks and valleys predicted by equations 4 and 

35 5 are observed in the experimental data; essentially, the intensity maxima are found 
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when the phase function is antisymmetric with' respect to the central wavelength of 
the pulse and minima when it is symmetric. ; 

[0081 ] More specifically, theoretical and experimental phase-mask control 
of two and three photon induced fluorescence is shown in Figures 14A - 14H. 
5 Equations 4 and 5 predict that as the phase mask is translated by an amount 8, the 
probability of two ( a P 2PA ") and three photon transitions ("P 3PA n ) is modulated, as 
illustrated in Figures 14A - 14D, for a half period sine mask (Figures 14A and 14B) 
and a full period sine mask (Figures 14C and 14D). Transform limited pulses yield a 
maximum value of 1. The small inserts in Figures 14A and 14C display the phase 

1 0 function at specific positions where maximum and minimum values of fluorescence 
take place (Figures 14E - 14H) wherein experimental two and three photon laser 
induced fluorescence measured for Coumarin and Stilbene, respectively, as a 
function of phase mask position 8 are shown. The phase masks used for these 
experiments were the same as those used in the calculations. Thus, the pulse 

15 shaping masks can be predetermined or fixed in shape based on calculations, 
experiments or learning program values for known equipment and known specimens. 
Example 2 

[0082] Experimental results for various samples obtained with a full- 
period sinusoidal phase mask are shown in Figures 15A - 15G. Figure 15A shows 
20 one photon laser induced fluorescence ("1 PLIF') of IR144 observed at 842 nm as a 
function of phase mask position. This measurement was made with 0.3 nJ/pulse to 
avoid nonlinear processes at the specimen. It is noteworthy that one photon process 
in the weak field regime show no dependence on phase shaping. Figure 15B shows 
results for the two photon laser induced fluorescence ("2PLIF") from Coumarin 
collected at 500 nm. The data in Figure 15C shows the dependence of 2PLIF in 
recombinant green fluorescent protein ("rGFP") detected at 505 nm. The data in 
Figure 15D corresponds to the intensity of the second harmonic generation ("SHG") 
signal at 405 nm from a 0.3mm p-barium borate crystal. The maximum and minimum 
signal for SHG coincides with that observed for 2PLIF but is not identical. 

[0083] With reference to Figure 15E, the dependence of three photon 
laser induced fluorescence ("3PLIF') from Trans-Stilbene is illustrated. The signal 
was collected at 350 nm as a function of 8. In this case, the maximum contrast 
(max:min) is measured to be 60:1 . The data in Figure 1 5F corresponds to the 3PLIF 
from Tryptophan residues in Con A, collected at 350 ( nm. In 3PLIF the maximum 
fluorescence signal is less than that obtained for transform limited pulses (when all 



25 
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the phases in the mask are set equal to zero), but the overall contrast ratio over the 
three-photon excitation is excellent, approaching two orders of magnitude. The data 
in Figure 15G corresponds to the continuum generation response (a nonlinear self- 
frequency modulation process yielding white light pulses) from a 3 mm slab of quartz 
5 detected at 600 nm. 

[0084] More specifically, Figures 15A - 15G demonstrate the 
experimental measurements of one and multi-photon emission obtained as a function 
of phase mask position 5. In all cases, the phase mask is a full period sine function. 
The signal measured with transform limited pulses is unity. The contrast ratio 

1 0 (max:min) is given in the upper right corner of each of the experimental plots. Here 
we find that the higher the order of the optical nonlinearity, the greater the contrast 
that we observe, therefore discrimination among different order processes is 
'possible. In other words, the higher the order, the greater the photons, which makes 
it easier for photon cancellation. Also, the greater the contrast ratio, the more the 

15 background noise is filtered out. 
Example 3 

[0085] Figure 16A presents the maximum discrimination between linear 
and nonlinear response observed for intense pulses (0.5 pJ/pulse). Separate 
detectors simultaneously collected the 1 PLIF from IRT44 solution and a portion of the 

20 continuum output. Maximum and minimum contrast ratios of >10 3 :1 and 1:0.6 were 
obtained for one photon process versus continuum, respectively, as shown in Figures 
16A and 16B. This control is extremely valuable when one is interested in linear 
processes under high-flux conditions, like in laser microscopy or in optical fiber 
communications. Using the simple phase function discussed earlier, particular 

25 windows of opportunity to control second versus higher order processes can be 
employed as demonstrated in Figures 16C and 16D. For certain values of 5, 
continuum generation even for relatively high intensity laser pulses (-1 pJ/pulse) can 
be completely suppressed. Figures 16C and 16D show that maximum and minimum 
contrast ratios of >10 3 :1 and 1:4 were obtained for 2PLIF versus continuum, 

30 respectively. 

[0086] Two photon transitions can be achieved while suppressing three 
photon processes for use in two photon microscopy or in two photon PDT. This type 
of control is much more difficult because once multiphoton transitions take place it is 
very difficult to stop at a particular order. A mixture of Coumarin and Fluoranthene 
35 were prepared to explore control of 2PLIF versus 3PLIF. Because fluorescence from 



W0 02/061799 PCT/US02/02548 



these two molecules overlaps the same spectral region, the separation between the 
two signals was achieved by temporal gating. Coumarin fluorescence was detected 
at 495 nm during the first 20 ns, while fluoranthene fluorescence was detected at 460 
nm with a gate that opened 40 ns after the initial rise and extended for 120 ns. 
5 Maximum and minimum contrast ratios of 1.4:1 and 1:2.2 were obtained for 2PLIF 
versus 3PLIF, respectively, as presented in Figures 16E and 16F. The contrast data 
presented in Figures 16A - 16F were obtained when transform limited pulses yielded 
equal intensities for the processes. Better contrast can be obtained using additional 
pulse shaping as described in the following section, especially as the multiphoton 

10 processes are detuned from resonance. 

Predetermined Pulse Shaping and Phase Control of Multiphoton Processes 

[0087] The present invention can utilize the presently invented" 
phenomenon of "Multiphoton Intrapulse Interference" as optimized for large 
molecules, proteins, and other condensed phase materials; through a combination of: 

15 (a) a chirped mask pulse shaper; and (b) a smooth function of phase versus 
frequency for the mask pulse shaper. The following formulas provide a predictive 
advantage for finding appropriate phase masks instead of using a learning program. 
The probability of two photon transitions can be calculated as follows for any given 
pulse shape: 

20 For an electric field with a carrier frequency o> 0 and a slow amplitude E 0 (t), 

oo 

[6] E(t) = Eo(t)e-toot and Eo(t) = -j=f Foi Q )e -iSlt dQ 

Y -oo 

where the Fourier image F 0 (G) around carrier frequency n=a>-a) 0 can be written as: 

oo 

^ F <>< Q ) = n=/Eo(t)e iQt dt, 

-OO 

the amplitude of two photon transition at resonance frequency co is: 
[8] A 2 (co) - jE(t) 2 e ia)t ilt = /EoWV^^tdt = / 4 (t) 2 e iAt dt) 

-°° -co ^ 

where detuning A=co-2o> 0 , the probability of two photon transition is: 
[9] P 2 (o» = |a 2 (g))| 2 . 
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[0088] Furthermore, the Fourier image of convolution is the product 
between Fourier images 
. HO] T(f*g) = (Tf)(Tg) 

where convolution (*, function from A) of two functions (f) and (g) is: 

[11] f*g=^7=/f(Q)g(A-Q)dQ. 

-oo 

Direct (T, function from Q) and inverse (T 1 , function from t) Fourier images are 
[12] T(f)=- 7 L/f(t)e i ^dtand T 1 (f) ^/f^e^da. 

5 Additionally, the relation between direct and reverse transforms is: 
[13] T 1 T(f) = TT 1 (f) = f. 

Thus, using the inverse transform, the formula can be written as: 
[1 4] f*g = r 1 T(f *g) = T 1 [(Tf)(Tg)] or 
formula [1 4] in integral form is as follows: 

00 00 00 00 

[1 5] Jf(£2)g(A-Q)d£2 = /e iAt [(-7=/f ^)e iQt d^(^fg(Q)e^dQ)]dt 

-oo -oo " -oo * -oo 

[0089] The time-frequency transformation can be calculated. Using the 
0 spectral presentation of formula [7] and convolution theorem of formula [15], formula 
[8] can be rewritten to obtain the formula for two photon transitions as follows: 

[1 6] A 2 (A) « /EoCtfe^dt - jF 0 (Q)F 0 (A-Q)dQ 
-°° -00 

This expression provides the two photon absorption amplitude given the spectrum of 
the laser pulse F 0 (Q) and the detuned spectrum of the F 0 (A-Q) that depends on the 
absorption spectrum of the sample. 

5 [0090] The probability of three photon transitions can be subsequently 

calculated. 

The complex amplitude of transition is: 



WO 02/061799 



PCT/US02/02548 



[1 7] A 3 (co) - /E(t) 3 e icot dt = jEoCtfe^dt, 

-CO _oo 

where detuning A=cd-3coo. Using the reverse Fourier presentation for the fields of 
formula [6], formula [17] can be rewritten as: 



co co 



[18] A 3 (©)°cpAtr/ Fo(Q)e 4Qt dQ ^^ 



Next, equation [18] can be rewritten using a new function G(Q) 



CO CO 



[1 9] A3(to) - Je' At [ jFo^e-^da, jG^e-^d^dt, 

-co -co -co 

where G(£2,) is defined as the kernel of the integral 

oo oa oo 

[20] fG(Q,)e-'^Q, = /Fo(Q 1 )e" i ^ t dn 1 Je^i^, and 

-oo -OO -co 

5 using the convolution formula [1 5], the following formula is obtained: 

co 

[21] A 3 (co)oc /FofpOGCA-QOdQL 

-CO 

[0091] The Fourier image of the Reverse Fourier image of equation [20] 
defines the intermediate function using relationship of equation [13] and the integral 
form of the convolution theorem expressed in formula [15] as: 

[22] "°° -°° -°° 

OO 

/Fo^FoCA-Qr^dQa 

-OO 

1 0 [0092] The final formula for the detuned A=co-3co 0 three photon transition 

is obtained by using equations [21] and [22] after changing the order of integration: 

oo 

oo OO 

[23] Aa(A) oc / Eo (t) 3 e^dt = / .../ F 0 (Q 1 )F 0 (Q 2 )F 0 (A-n 1 -Q 2 )dQ 1 dQ 2 

-OO 

such that the probability is: 



-oo -oo 
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[24] P 3 ((0) = |A 3 (0))| 2 . 

The method described above gave the formula for the n-photon transition by 
recurrence: 



oo oo 



[25] A„(A) oc jEo(t) n e iAt dt « J ... j F 0 (Q,). . .Fo^OFoOa . . .-Q^dQ, . . .dQo., 

where detuning is A=co-nco 0 . Thus, 
[26] P n (o>) oc |A n (co)| 2 . 

5 [0093] It is also desirable to take into account inhomogeneous 

broadening (as encountered in solutions and condensed phase materials). The 
integrated probability for the n-photon transition in molecules with spectral a density 
g n (o)) with amplitude defined by formula [25] is proportional to the weighed average 

oo 

[27] P»= /&(0))|A n (CD)| 2 d0). 

— OO 

Normalization for the case of transform limited laser pulse is N n and 

oo 

[28] N,= / 6 (a))|A llI (a))pdcD b - 

— OO 

10 where 

OO 

[29] A TXn (to)= . . . / Pb(Q,)|. . . jFoCA-OHFoCA-Q! . . .-G^dO! . . .dQ^. 



[0094] The preceding formulas [6] - [29] give the general result. The 
following parameters must be defined, however, for a user to define a phase mask 
that would minimize or maximize a particular multiphoton process. First, the laser 

15 pulse spectrum of Figure 18A must be defined. The shorter the pulses (broader 
spectrum), the better the control. 45 fs pulses have been satisfactorily used but 20 or 
10 fs would lead to even better results. The carrier frequency (or center wavelength) 
must also be defined by availability. Tuning the wavelength of the pulse could 
enhance certain processes but is not typically required. Secondly, the phase 

20 modulator (or alternately, the SLM, deformable mirror, chirped mirror, etc.) should 
cover the entire pulse spectrum and must be defined. Thirdlv. a Dhase mask 
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definition should be introduced. The simple sine function of Figure 18B works 
remarkably well, yet other functions that can become symmetric and antisymmetric 
as a function of their position are also suitable. Fourthly, the addition of positive or 
negative linear chirp p further enhances the observed control, as expressed in Figure 
18C, and should be defined. The phase mask used in the examples presented 
herein is defined by 



t 3 °3 q> m (U) = <Min 



J 



where 8 is the position of the sine function (centering) across the spectrum, <|> a is the 
maximum phase delay, and Np ixe i is the number of pixels in the SLM, as illustrated in 
Figure 18B. 

10 [0096] When chirp is added, it can be defined by 

[31] <k(Q) = y 2 pQ 2 

where p is the amount of linear chirp expressed in Figure 18C. Thus, the complete 

phase mask with chirp is: 

[32] M) = UU)-H>c(X). 

[0097] Figures 19A and 19B shpw the calculated two and three photon 
absorption probability using the equation presented and the absorption spectra as 

15 calculated by the dotted lines in Figures 20A - 20C. Figure 20C shows the 
calculated ration two:three photon absorption for the two different combinations of 
absorption spectra given in Figures 20A and 20B. Accordingly, robust control of 
multiphoton processes in molecules, proteins and nonlinear optical materials can be 
achieved through either adaptive, active and self optimizing learning programs and 

20 control systems, or through calculated, predetermined or fixed, passive laser beam 
pulse shaping devices. Therefore, inexpensive fixed phase masks can be designed 
before the experiment, and even without computer controlled shapers and learning 
programs, to control the order of multiphoton processes for large, complex 
molecules, proteins in photodynamic therapy, optical tomography, surgery (such as 

25 laser cutting by five or greater photon wave conveyance to maximize nonlinear 
energy), and photochemistry control of, for example: (a) photopolymerization (by 
photon pair switching to seed the process), (b) charge transfer, (c) radical reaction, 
(d) nucleophelic attack and (e) electrophylic attack. 
Communications 

30 [0098] With reference to Figure 21, a fifth preferred embodiment of a 

la.Q^r OvHtatinn onnorotnn Ar\A .x iu- * : | , 
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423, an optical fiber 451, a laser beam pulse shaper device 429, a laser beam pulse 
un-shaper device 453, and a receiver 441 which includes an optical switch or sensor 
and the related circuitry and electrical control unit. Laser 423 emits a series of laser 
beam pulses, each shorter than 1 ps, into the connected fiber 451 . Pulse shaper 
5 device 429 is of a predetermined mask type with a fixed pulse characteristic varying 
shape (such as with calculated sine wave surface shapes) and has three elements 
connected to fiber 451: a dispersive element 455 such as a fiber that incorporates a 
diffraction grating; a phase mask element 457 that can be made using a doped glass 
or polymer sheet; and a dispersive element 459, like element 455 but reversed, for 

1 0 accepting spectrally dispersed light and coupling it back to fiber 451 . 

[0099] The shaped laser beam pulse is capable of traveling long 
distances through fiber 451 without suffering nonlinear distortion because of the 
unique phase function imprinted or formed on shaper device 429. For example, the 
red color spectrum may be advanced in front of the blue color spectrum in a precise 

15 sine manner. Un-shaper device 453 subsequently reverses the phase changes 
introduced by shaper device 429. It is constructed the same as the shaper device 
but with a different phase mask element 461 that compensates for the pulse 
characteristic changes made by mask element 457. Alternately, an acousto-optic 
modulator or transient grating can be used for optical switching through constructive 

20 or destructive reference of waves. Shaping and unshaping can also be 
accomplished by means of a chirped mirror or spectral masks. 

[00100] Thus, the present invention's ability to precisely control the laser 
beam pulse shape or other characteristic, especially for nonlinear or multiphoton 
emissions, significantly improves the quality of the communication transmission while 

25 minimizing self-focusing, self phase modulation and possible destruction of the fiber. 
The pulse characteristic control of ultrafast laser beam pulses, as described in all of 
the embodiments herein, should minimize, if not prevent, multiplicative noise effect 
disruption of nonlinear propagation channels in fiberoptic lines, as discussed in Mitra, 
et at., "Nonlinear Limits to the Information Capacity of Optical Fibre 

30 Communications," Nature, vol. 411, pp. 1027-1030 (June 28, 2001). It is further 
envisioned that this type of pulse shaping apparatus can be employed within salt 
water oceans for submarine-to-submarine communications using shorter than 1 ps 
laser pulses. 

[00101] While the preferred embodiment of the control system and 
35 apparatus of the present invention have been disclosed, it should be appreciated that 
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various modifications can be made without departing from the spirit of the present 
invention. For example, other laser beam pulse characteristics can be varied and 
employed with the present invention beyond the pulse shaping and duration 
characteristics preferably described. Furthermore, additional software subroutines 
5 and statistical analyses can be employed. Moreover, other optical and pulse shaping 
devices can be used in place of those described, such as deformable mirrors and the 
like. Finally, analog, solid state and fiber optic electrical control circuits can be 
substituted for or used in addition to a microprocessor and other computer circuitry. 
Various optics, including lenses and mirrors, can be used to. achieve collimation or 

10 focusing. Additionally, dispersive optics, such as gratings and prisms, can be 
interchanged. Detection of the laser induced processes may use various 
spectroscopic methods including laser induced fluorescence, Raman spectroscopy, 
nuclear magnetic resonance, gas chromatography, mass spectrometry and 
absorbtion spectroscopy. While various materials, specimens and components have 

15 been disclosed, it should be appreciated that various other materials, specimens, and 
components can be employed. It is intended by the following claims to cover these 
and any other departures from the disclosed embodiments which fall within the true 
spirit of this invention. : 
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The invention claimed is: 

1. An apparatus for use with a specimen, the apparatus 

comprising: 

5 a femtosecond laser operable to emit a laser beam of less than about 

50 femtosecond pulse duration upon the specimen; 

a pulse shaper operable to shape the laser beam; 

a device operable to detect characteristics of the specimen when the 
shaped laser beam has been emitted upon the specimen; and 
10 an electrical control system connected to the device and the pulse 

shaper, the electrical control system operably varying the pulse shaping performance 
of the pulse shaper for subsequent laser beam emissions based on an evaluation of 
the characteristics detected by the device. 

2. The apparatus of Claim 1 wherein the pulse shaper includes a 

15 Fourier plane. 

3. The apparatus of Claim 2 wherein the pulse shaper further 
includes at least one collimating optic and a dispersive optic. 

4. The apparatus of Claim 3 wherein the pulse shaper includes a 
second collimating optic and a second dispersive optic, wherein the Fourier plane is 

20 optically located between the collimating optics. 

5. The apparatus of Claim 1 wherein the specimen is a protein. 

6. The apparatus of Claim 5 wherein the electrical control system 
determines the protein sequence of the specimen. 

7. The apparatus of Claim 5 wherein the laser beam selectively 
25 cleaves bonds in the specimen. 

8. The apparatus of Claim 5 wherein there is direct laser 
desorption of the protein specimen free of a matrix. 

9. The apparatus of Claim 1 wherein the device is a mass 

spectrometer. 

30 10- The apparatus of Claim 9 wherein matrix-assisted laser 

desorption ionization, post-source decay mass spectrometry is employed to analyze 
the specimen by the mass spectrometer and electrical control system. 

11. The apparatus of Claim 1 wherein the pulse shaper includes 
an acousto-optic modulator operable to control amplitude and phase of the laser 

35 beam. 
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12. The apparatus of Claim 11 wherein the acousto-optic 
modulator includes an anti-reflection coated Tellurium Dioxide crystal with a piezo- 
electric transducer. 

13. The apparatus of Claim 1 wherein the pulse shaper includes a 
5 liquid crystal display operable to control amplitude and phase of the laser beam. 

14. The apparatus of Claim 1 wherein the pulse shaper includes a 
deformable mirror. 

15. The apparatus of Claim 1 wherein the electrical control system 
automatically determines the most desirable laser beam pulse shape based on prior 

1 0 device results using substantially randomly employed pulse shapes. 

16. The apparatus of Claim 1 wherein the femtosecond laser 
operably transmits a laser beam of less than 1 1 femtosecond duration. 

17. The apparatus of Claim 1 wheVein the specimen is multi- 

15 molecular. 

18. A matrix-assisted laser desorption ionization system 
comprising: 

a multi-molecular specimen; 

a laser operable to emit a laser beam upon the specimen; 
20 a Pu'se shaper operable to shape the laser beam; 

a mass spectrometer operable to detect characteristics of the 
specimen if the shaped laser beam has been emitted upon the specimen; and 

an electrical control system connected to the mass spectrometer and 
the pulse shaper, the electrical control system automatically varying the pulse shapes 
25 for subsequent laser beam emissions based on at least one signal generated by the 
mass spectrometer. 

1 9. The system of Claim 1 8 wherein the specimen is a protein. 

20. The system of Claim 19 wherein the electrical control system 
determines the protein sequence of the specimen. 

60 21 • The s y ste m of Claim 19 wherein the laser beam selectively 

cleaves bonds in the specimen. 

22. The system of Claim 19 wherein there is direct laser 
desorption of the protein specimen free of a matrix. 

23. The system of Claim 18 wherein the pulse shaper includes at 
35 least one collimating optic and a grating. 
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24. The system of Claim 23 wherein the pulse shaper includes a 
second optic and a second grating. 

25. The system of Claim 18 wherein the pulse shaper includes an 
acousto-optic modulator operable to control amplitude and phase of the laser beam. 

5 26. The system of Claim 25 wherein the acousto-bptic modulator 

includes an anti-reflection coated Tellurium Dioxide crystal with a piezo-electric 
transducer. 

27. The system of Claim 18 wherein the laser is a femtosecond 
laser operable to emit a laser beam of less than . about 21 femtosecond pulse 

1 0 duration upon the specimen. 

28. The system of Claim 18 wherein the pulse shaping is achieved 
by a set of interchangeable passive phase amplitude masks. 

29. The system of Claim 18 wherein matrix-assisted laser 
desorption ionization, post-source decay mass spectrometry is employed to analyze 

1 5 the specimen by the mass spectrometer and electrical control system. 

30. The system of Claim 18 wherein the pulse shaper includes a 
liquid crystal display operable to control amplitude and phase of the laser beam. 

31 . The system of Claim 1 8 wherein the pulse shaper includes a 
deformable mirror. 

20 32. The system of Claim 18 wherein the electrical control system 

automatically determines the most desirable laser beam pulse shape based on prior 
device results using substantially randomly employed pulse shapes. 

33. The system of Claim 18 wherein the laser is a femtosecond 
laser which operably transmits a laser beam of less than 1 1 femtosecond duration. 
25 34. The system of Claim 18 wherein: 

the electrical control system operably causes the laser and pulse 
shaper to transmit an initial random population of pulse shapes upon the specimen; 

the electrical control system operably compares a target to each pulse 
shape result as sensed by the mass spectrometer; 
30 the electrical control system then operably causes the laser to emit 

laser beam pulses and the pulse shaper to transmit subsequent sets of varying pulse 
shapes the results of which are then compared to at least one of the pulse shape 
results of the prior iterations; and 

the electrical control system operably determines if statistical 
35 convergence has been obtained from repeated pulse shape set iterations. 
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55. The system of Claim 54 wherein matrix-assisted laser 
desorption ionization, post-source decay mass spectrometry is employed to analyze 
the specimen by the mass spectrometer and electrical control system. 

56. The system of Claim 46 wherein the specimen is a protein. 

5 57. The system of Claim 46 wherein the electrical control system 

automatically determines the most desirable laser beam pulse shape based on prior 
device results using substantially randomly employed pulse shapes. 

58. The system of Claim 46 wherein the femtosecond laser 
operably transmits a laser beam of less than 1 1 femtosecond duration. 
1 0 59. A system for use with a sample, the system comprising: 

laser beam pulses operably emitted upon the sample; and 
a pulse shaper operable to vary shapes of the laser beam pulses; 
wherein the shaped laser beam pulses chemically modify the sample. 

60. The system of Claim 59 wherein the shaped laser beam pulses 
1 5 selectively cleave atomic bonds in the sample. 

61. The system of Claim 60 wherein chemical bonds are 
selectively cleaved in the sample. 

62. The system of Claim 59 further comprising a laser operably 
emitting the laser beam pulses with less than 50 femtosecond durations. 

20 63 - The system of Claim 59 further comprising an electrical control 

system operably causing the pulse shaper to vary the laser beam pulse shapes 
based on modification results in an automatic manner. 
64. A system comprising: 
a specimen substrate having isolated molecules; 
25 at least one laser beam pulse operably emitted upon the molecules; and 

a pulse modification device operably changing an excitation characteristic of the 
emitted pulse; 

wherein the 'modified laser beam pulse causes desorption of the 
isolated molecules. 

30 65 - The system of Claim 64 wherein the specimen does not use a 

matrix between the molecules, and the pulse modification device is a pulse shaper. 
66. A system for use with a sample, the system comprising: 
a laser operable to emit a laser beam pulse upon the sample; 
a mass spectrometer operable to analyze ionization of the sample- 

35 and 



WO 02/061799 



PCT/US02/02548 



an electrical control system operable to cause the laser beam pulse to 
be modified prior to it being received by the sample, the electrical control system 
further modifying and substantially optimizing at least one characteristic of 
subsequent laser beam pulses. 
5 67. The system of Claim 66 wherein the modification and 

optimization includes a laser beam pulse having a combination of pulse durations 
determined by the electrical control system. 

68. The system of Claim 66 wherein the modification and 
optimization includes a laser beam pulse having a combination of pulse shapes 

10 determined by the electrical control system. 

69. The system of Claim 66 wherein the modification and 
optimization includes a laser beam pulse having a combination of pulse wavelengths 
determined by the electrical control system. 1 

70. A system comprising: 
1 5 a protein specimen; 

a laser beam pulse of less than about 50 femtosecond duration 

operably emitted upon the specimen; 

a device operable to detect characteristics of the specimen when the 

laser beam is emitted upon the specimen; and 
20 an electrical control system connected to the device, the electrical 

control system being operable to sequence proteins for the specimen based at least 
in part from signals generated by the device. 

71 . The system of Claim 70 further comprising a pulse shaper 
operable to vary at least one of the amplitude and phase of the laser beam. 

25 72. The system of Claim 71 wherein the pulse shaper includes a 

Fourier plane, at least one collimating mirror and a grating. 

73. The system of Claim 71 wherein the pulse shaper uses a 
smooth phase function across the spectrum. 

74. The system of Claim 70 further comprising an acousto-optic 
30 modulator operable to control amplitude and phase of the laser beam. 

75. The system of Claim 70 further comprising a liquid crystal 
display operable to control amplitude and phase of the laser beam. 

76. The system of Claim 70 further comprising a deformable mirror 

operably to assist in shaping the laser beam pulse. 
35 77. The system of Claim 70 wherein the specimen is a protein. 
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35. An apparatus for use with a specimen, the apparatus 

comprising: 

a laser operable to emit a laser beam; 
a pulse shaper operable to shape the laser beam; 
a mass spectrometer operable to sense characteristics of the 
specimen if the shaped laser beam has been emitted upon the specimen; and 

an electrical system connected to the. mass spectrometer and the 
pulse shaper, the electrical system operably varying the later emitted laser beam 
pulse shapes based on a comparison of the mass spectrometer results, the electrical 
system operably causing the laser to emit at least 1 ,000 laser pulses at the specimen 
and the specimen providing desirable ionization substantially throughout all of the 
pulses without requiring a replacement of the specimen. 

36. The apparatus of Claim 35 wherein the specimen is a protein. 

37. The apparatus of Claim 36 wherein the electrical system 
15 determines the protein sequence of the specimen. 

38. The apparatus of Claim 36 wherein the laser beam selectively 
cleaves bonds in the specimen. 

39. The apparatus of Claim 36 wherein there is direct laser 
desorption of the protein specimen free of a matrix. 

20 40 - The apparatus of Claim 35 wherein matrix-assisted laser 

desorption ionization, post-source decay mass spectrometry is employed to analyze 
the specimen by the mass spectrometer and electrical system. 

41. The apparatus of Claim 35 wherein the electrical system 
automatically determines the most desirable laser beam pulse shape based on prior 
mass spectrometer results using substantially randomly employed pulse shapes. 

42. The apparatus of Claim 35 wherein the laser is a femtosecond 
laser operably providing laser beam pulses of less than 21 femtosecond duration. 

43. The apparatus of Claim 35 wherein the laser is a femtosecond 
laser which operably transmits a laser beam of less than 1 1 femtosecond duration. 

3t 44 ' The a PParatus of Claim 35 wherein the specimen is multi- 

molecular. 

45. The apparatus of Claim 35 wherein: 
the electrical system operably causes the laser and pulse shaper to 
transm.t an initial random population of pulse shapes upon the specimen; 



25 
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the electrical system operably compares a target to each pulse shape 
result as sensed by the mass spectrometer; 

the electrical system then operably causes the laser to emit laser 
beam pulses and the pulse shaper to transmit subsequent sets of varying pulse 
5 shapes the results of which are then compared to at least one of the pulse shape 
results of the prior iterations; and 

the electrical system operably determines if statistical convergence 
. has been obtained from repeated pulse shape set iterations. 
46. A system for conducting a test, the system comprising: 
10 a multi-molecular specimen; 

a laser beam operably emitted upon the specimen; 
a device operable to detect excitation characteristics of the specimen 
when the laser beam is emitted upon the specimen; and 

an electrical control system connected to the device, the electrical 
1 5 control system being operable to vary a feature of the emitted laser beam to cleave 
an atomic bond in the specimen. 

47. The system of Claim 46 further comprising a pulse shaper 
operable to vary at least one of the amplitude and phase of the laser beam. 

48. The system of Claim 47 wherein the pulse shaper further 
20 includes at least one collimating optic and a dispersive optic. 

49. The system of Claim 48 wherein the pulse shaper includes a 
second collimating optic and a second dispersive optic, and a Fourier plane is 
optically located between the collimating optics. 

50. The system of Claim 46 wherein pulse shaping is achieved 
25 with interchangeable phase amplitude altering optics. 

51 . The system of Claim 47 wherein the pulse shaper includes an 
acousto-optic modulator operable to control amplitude and phase of the laser beam. 

52. The system of Claim 47 wherein the pulse shaper includes a 
liquid crystal display operable to control amplitude and phase of the laser beam. 

30 53. The system of Claim 47 wherein the pulse shaper includes a 

deformable mirror. 

54. The system of Claim 46 wherein the device is a mass 
spectrometer and the identify of the specimen was unknown prior to the test. 
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78. The system of Claim 70 wherein the device is a mass 

spectrometer. 

79. The 'system of Claim 70 wherein matrix-assisted laser 
. desorption ionization, post-source decay mass spectrometry is employed to analyze 

5 the specimen. 

80. The system of Claim 70 wherein the electrical control system 
automatically determines the most desirable laser beam pulse shape based on prior 
device results using substantially randomly employed pulse shapes. 

81. A control system for use with a specimen, the system 

10 comprising: 

a laser operable to emit a laser beam upon the specimen; 
a pulse shaper operable to shape the laser beam; 
a detection device operable to detect characteristics of the specimen 
when the shaped laser beam has been emitted upon the specimen; and 

15 a control device connected to the detection device and the pulse 

shaper; 

the control device automatically, varying the puise shapes for 
subsequent laser beam emissions based at least in part on at least one signal 
generated by the detection device; f 
20 tne control device operably comparing a target to each pulse shape 

result as detected by the detection device; 

the control device operably causing the laser to emit laser beam 
pulses and the pulse shaper.to transmit subsequent sets of varying pulse shapes, the 
results of which are then compared to at least one of the pulse shape results of the 
25 prior iterations; and 

the control device operably determining if statistical convergence has 
been obtained from repeated pulse shape set iterations. 

82. The system of Claim 81 wherein the laser operably transmits a 
laser beam of less than about 20 femtosecond duration. 

30 83 - The system of Claim 81 wherein the detection device is part of 

a matrix-assisted laser desorption ionization, post-source decay system and the 
detection device is a mass spectrometer, and the control device includes a 
microprocessor. 

84. The system of Claim 81 wherein the specimen is a protein and 
35 the control device assists ^determining a sequence of the protein. 
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85. The system of Claim 81 wherein the control device varies the 
laser beam emission in an adaptive manner. 

86. A system for use with living tissue, the system comprising: 
a high peak intensity laser beam pulse; and 

5 a device operable to change a characteristic of the pulse prior to 

. emission of the pulse upon the living tissue; 

wherein nonlinear transitions induced by each pulse are controlled. 

87. The system of Claim 86 wherein the device uses a pulse 
shaper and the desired eradicated substances in the tissue undergo two photon 

10 absorbtion. 

88. The system of Claim 86 wherein the pulse has a duration of 
less than fifty one femtoseconds. 

89. The system of Claim 86 further comprising generating an 
optical tomography image produced by the shaped pulse passing through the tissue. 

15 90. The system of Claim 86 wherein the pulse shaper enhances 

two photon absorbtion by a therapeutic substance and substantially prevents three 

photon induced damage of adjacent healthy tissue. 

91 . The system of Claim 86 wherein the device includes a chirped 

phase mask modifying the beam. 
20 . 92. The system of Claim 86 wherein the pulse is shaped to 

enhance targeted mutliphoton damage to modify or destroy certain molecules in the 

living tissue. 

93. A communications system comprising: 
a laser operably emitting laser beam pulses; 

25 a first device operable to modify the phase of the pulses emitted from 

the laser; 

and 

30 distortion. 

94. The communications system of Claim 93 further comprising at 
.. least one of the following transmissive media: (a) an optical fiber; (b) air; and (c) 

water; connecting the first and second devices. 

95. The communications system of Claim 93 wherein the first and 
35 second devices are laser beam pulse shapers. 



a second device operable to substantially reverse the modified pulses; 
a receiver operably receiving the reversed pulses substantially free of 
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CONTROL SYSTEM AND APPARATUS FOR 
USE WITH LASER EXCITATION OR IONIZATION 

BACKGROUND AND SUMMARY OF THE INVENTION 
5 [0001] The present invention generally relates to control systems and 

apparatuses for use with laser excitation or ionization, and more particularly to a 
system and apparatus which employs a laser, pulse shaper, mass spectrometer and 
electrical control system. 

[0002] Conventionally, laser desorption mass spectrometry has been 
10 used with a fixed laser beam pulse shape and computers for simple chemical 
analysis processes on purified molecules with or without a matrix. The laser beam 
pulse shape was not considered an important parameter and was not modified; 
whatever fixed shape was set by the manufacturer for the ultraviolet laser was used 
in the tests. The general concept of typically laser selective ion formation from 
15 molecules in a molecular b^am is disclosed in the following publication: Assion et al., 
"Control of Chemical '. Reactions by Feedback-Optimized Phase-Shaped 
Femtosecond Laser Pulses," Science, Vol. 282, page 919 (October 30, 1998). The 
pulse shaping process with a learning algorithm is disclosed in Judson et al., 
Teaching Lasers to Control Molecules," Physical Review Letters, Vol. 68, No. 10, 
20 page 1500 (March 9, 1992). It is noteworthy, however, that the Assion article 
discloses use of an 80 femtosecond laser pulse and requires molecules to be 
isolated in a molecular beam, while the Judson article discloses use of a one 
nanosecond laser pulse and is purely conceptual as it does not include experimental 
results. 

25 [0003] It is also known to employ nanosecond lasers for matrix-assisted 

laser desorption ionization (hereinafter "MALDI"). Examples of this are disclosed in 
U.S. Patent No. 6,130,426 entitled "Kinetic Energy Focusing for Pulsed Ion 
Desorption Mass Spectrometry" which issued to Laukien et al. on October 10, 2000, 
and U.S. Patent No. 6,111,251 entitled "Method and Apparatus for MALDI Analysis" 

30 which issued to Hillenkamp on August 29, 2000; both of these patents are 
incorporated by reference herein 

[0004] Until recently, commercially practical femtosecond lasers have 
been unavailable. For example, lasers which can generate 10 femtosecond or less 
laser pulse durations have traditionally been extremely expensive, required 

35 unrealistically high electrical energy consumption (for extensive cooling, by way of 
example) and denendeH 
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thereby leading to commercial impracticality. The efficiency of sub-10 femtosecond 
lasers was not practical until the year 2000 because of the prior need for dyes and 
flash lamps instead of YAG and Ti: Sapphire crystals pumped by light or laser 
emitting diodes. 

5 [0005] Furthermore, the traditional role of the laser in a mass 

spectrometer with MALDI is to provide energy to the matrix molecules, wherein this 
energy dissipates and causes evaporation and ionization of the protein analyte 
dissolved in it. The laser, therefore, plays an indirect role that depends on energy 
transfer processes that may take from picoseconds to microseconds. Because 

10 excitation is indirect, pulse wavelength has not been found to cause significant 
differences in the outcome. Direct laser excitation of the proteins with nanosecond 
. lasers typically causes the proteins to char. 

[0006] In contrast, the present invention uses a different approach to 
MALDI in which the laser plays a more active and direct role in the ionization and 

15 even selective fragmentation of the analyte proteins.^; Shaped femtosecond pulses 
are required to achieve this goal. The optimum pulse shape cannot be found using 
the traditional laser sources, and trial and error. This is because the search for an 
optimal laser pulse shape involves a very wide range of possibilities. For example, if 
a 1 00 femtosecond laser pulse is used to produce pulse trains as long as several 

20 picoseconds in duration, splitting the laser beam spectrum into at least 100 spectral 
components is required since the length of the pulse is roughly inversely proportional 
to the band width. Since each component can be attenuated in 10 steps or phases 
shifted over 10 angles, then there are (10 X10) 100 different possible pulse shapes, 
and it would be impractical to systemically explore even a subset of these pulse 

25 shapes through conventional trial and error methods. 

[0007] Laser induced, selective chemical bond cleavage has also been 
explored but with fairly limited success. It is believed that very simple molecules, 
such a HOD (partially deuterated water), have had only the OH and OD bonds 
cleaved with a nanosecond narrow line laser to vibrational^ excite the specimen and 

30 tf n an ultraviolet laser pulse was employed to perform the cleaving. The desired 
laser frequency for vibrational excitation could be determined a priori in the gas- 
phase sample. More importantly, the HOD molecule is unique because the energy 
can be deposited in one of the bonds and it remains there for very long times, which 
are longer than nanoseconds. For the HOD experiments using selective bond 

35 excitation, no appreciable pulse shaping was used. This method was not known to 
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have been employed for a protein or MALDI process, and was not known to have 
been successfully used for any other atomic bonds in other molecules, especially not 
in a condensed phase. It is also noteworthy that MALDI, with a matrix, has been 
used in an attempt to perform limited bond cleavage, as is discussed in U.S. Patent 
5 No. 6,156,527 entitled "Characterizing Polypeptides" which issued to Schmidt et al. 
on December 5, 2000, and is incorporated by reference herein. However, the 
approach of Schmidt et al. does not modify and optimize the laser pulse shape or 
other laser properties to achieve iimited bond cleavage. 

[0008] In accordance with the present invention, a control system and 

1 0 apparatus for use with laser excitation or ionization is provided. In another aspect of 
the present invention, the apparatus includes a laser, pulse shaper, detection device 
and control system. A further aspect of the present invention employs a 
femtosecond laser and a mass spectrometer. In yet another aspect of the present 
invention, the control system and apparatus are used in a MALDI process. Still 

1 5 another aspect of the present invention employs the control system and apparatus to 
cleave chemical bonds in a specimen and/or to determine the amino acid sequence 
of a protein specimen. Photodynamic therapy and fiber optic communication 
systems use the laser excitation apparatus with additional aspects of the present 
invention. A method of ionizing and determining a characteristic of a specimen is 

20 also provided. 

[0009] The control system and apparatus of the present invention are 
advantageous over conventional constructions since the present invention allows for 
analysis and identification of constituents of complex and unknown molecules, such 
as those used in a MALDI process or proteins, in a relatively quick and automated 

25 manner. The present invention advantageously determines optimum laser conditions 
for maximizing the sensitivity of MALDI based protein sequencing, and to examine 
ion formation efficiencies for various matrices using tailored laser pulses. The 
present invention is also advantageously used to control the degree and type of . 
fragmentation for automated protein sequencing. Furthermore, the adaptive laser 

30 source permits the optimal desorption from an insoluble protein source and allows for 
ionization analysis of a protein with or without a matrix. 

[0010] The present invention is advantageous by employing ultra-fast 
laser beam pulses which can be repeatedly transmitted onto a specimen at least 
1,000 times without replacing the specimen and without significant degradation of 

35 results. The ultra-fast laser also does not over-heat or "cook" a specimen, such as a 
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protein. Recent improvements and efficiencies of femtosecond lasers have allowed 
for their commercially practical usefulness with the present invention. The automated 
feedback and pulse shaping of the control system of the present invention enhances 
signal-to-background sensitivity, especially for MALDI-based protein sequencing, 
5 while also statistically optimizing the process; this leads to significant time, cost and 
accuracy improvements. Additional advantages and features of the present invention 
will become apparent from the following description and appended claims, taken in 
conjunction with the accompanying drawings. 

10 BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] Figures 1 through 3. are diagrammatic views showing a first 
preferred embodiment of a control system and apparatus of the present invention; 

[001 2] Figure 4 is a flow chart showing the operation of the first preferred 
embodiment control system and apparatus; 
15 [0013] Figure 5A through 7B are sets of laser beam pulse shapes 

employed with the first preferred embodiment control system and apparatus; 

[0014] Figure 8 is an exemplary ionization fragmentation chart employed 
with the first preferred embodiment control system and apparatus; 

[0015] Figures 9 through 10B are flow charts for the method and 
20 computer software operation employed with the first preferred embodiment control 
system and apparatus; 

[0016] Figure 11 is an exemplary molecular structure cleaved by the first 
preferred embodiment control system and apparatus; 

[0017] Figure 12 is a diagrammatic view showing a second preferred 
25 embodiment of a laser excitation apparatus of the present invention; 

[0018] Figures 13A-13C are schematic and graphical representations of 
two photon and three photon induced fluorescence employed with the laser excitation 
apparatus; 

[0019] Figures 14A-14H are sets of laser beam pulse shapes employed 
30 with the laser excitation apparatus for two and three photon induced fluorescence; 

[0020] Figures 15A-15G are sets of laser beam pulse shapes employed 
with the laser excitation apparatus for two and three photon induced fluorescence; 

[0021] Figures 16A-16F are sets of pie charts and laser beam pulse 
shape graphs showing contrast ratios of the laser excitation apparatus; 
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[0022] Figure 17 is a diagrammatic view showing third and fourth 
preferred embodiments of the present invention laser excitation apparatus applied to 
optical coherence tomography and photodynamic therapy; 

[0023] Figures 18A - 18C are graphs showing the laser beam pulse 
5 spectrum employed with the laser excitation apparatus; 

[0024] Figures ^9A and 19B are graphs showing the calculated two and 
three photon absorption probability employed with the laser excitation apparatus; 

[0025] Figures 20A - 20C are graphs showing the calculated two and 
three photon absorption probability employed with the laser excitation apparatus; and 
10 [0026] Figure 21 is a diagrammatic view showing a fifth preferred 

embodiment of the present invention laser excitation apparatus applied to fiber optic 
communications. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

1 5 [0027] The first preferred embodiment of a control system and apparatus 

21 of the present invention for use with laser excitation or ionization is generally 
shown in Figures 1 and 2. Apparatus 21 includes a femtosecond laser 23, an 
upstream grating 25, an upstream convex mirror 27, a laser beam pulse shaper 29, a 
downstream concave mirror 31, a downstream grating 33, a matrix-assisted laser 

20 desorption ionization device 35, and a personal computer 37. Personal computer 37 
has a microprocessor based electrical control system, an output screen, a data 
storage device, an input keyboard, and a removable disk. More specifically, the 
MALDI device provides a time-of-flight mass spectrometer ('TOF MS") 39. A sample 
or specimen 41 to be analyzed is placed within mass spectrometer 39. Bursts or 

25 pulses of a laser beam 43 are emitted from laser 23, through the optics 25, 27, 31 
and 33, as well as through pulse shaper 29, and onto specimen 41; this causes 
fragmentation and ionization of a top layer of the specimen for detection and sensing 
by mass spectrometer 39 for further evaluation, analysis, comparison and 
subsequent control by personal computer 37. 

30 [0028] The laser is preferably an ultra-fast femtosecond and high peak 

intensity (with a typical peak greater than 10 10 watts/cm 2 ) laser which preferably 
emits laser beam pulses of less than 100 femtosecond duration, and more preferably 
at or less than .50 femtoseconds, and for certain applications (such as, but not limited 
to, sequencing) even more preferably at or less than 10 femtosecond duration, for 

35 each pulse burst or shot. The intense optical pulses that are required to modifv 
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material are formed in a Kerr-Lens modelocked titanium sapphire oscillator. Such 
lasers are capable of producing hundreds of nanometers of coherent 'bandwidth, 
although only about 10 nm are typically used. The output is amplified in a 1 kHz 
regenerative chirped pulsed amplifier. The output pulse is typically 100 fs .long with a 
5 central wavelength of out 790 nm and total pulse energy of 0.1 to 1 mJ. Preferred 
lasers include: the Hurricane model from Spectra Physics Inc., which is diode 
pumped and gives 0.8 mJ per pulse with sub-50 fs pulses at 1 kHz; and the CPA- 
2001+ model from Clark-MXR Inc., which gives 1.3 mJ per pulse with sub-150 fs 
pulses at 1 kHz, pumping a Clark-MXR Inc. non-collinear parametric amplifier 
10 (hereinafter "NOPA") which produces 0.2 mJ per pulse, and is capable of generating 
sub-20 fs pulses. This NOPA system can even produce pulses between 10 fs and 
4.5 fs. 

[0029] A VESTEC 2000 MALDI TOF mass spectrometer is believed to be 
suitable for this invention, although most commercial MALDI instruments can be 

15 adapted with the femtosecond laser, pulse shaper and feedback learning control 
method described herein. During extraction, all of the ions obtain the same energy in 
the 30-kV ion acceleration region, and because K.E. = Vz mv 2 , the lightest ions 
achieve the highest velocity and, thus, reach the detector first. This transient (for 
example, having a duration of 300 microseconds) mass spectrum is recorded by a 

20 transient recorder at the detector. It is common practice to sum many (10-100) of 
this transient mass spectra to produce a sound spectrum from an ion-counting 
statistics criterion. 

[0030] In Figure. 3, a conceptual schematic is shown outlining the four 
probable mechanisms theorizing where the proton transfer takes place. The laser 

25 impinges on matrix crystal specimen 41 creating a plume of highly charged particles. 
Some of the particle are proton rich, such as M2H + and are good proton donors. Gas 
phase collisions with the large and slow moving proteins may lead to the formation of 
the protonated analyte. Other gas-borne particles include clusters containing the 
naturally charged analyte and matrix molecules. As the clusters evaporate, the 

30 charged analyte remains. Alternatively, the laser causes local melting of the crystal. 
The proton transfer occurs in the liquid phase between the excited state (high acidity) 
matrix molecules and the analyte. Finally, the proton transfer could also occur in the 
solid phase. The. excitons created by the laser near the protein can relax by reactive 
processes such as proton transfer. At present, it is unknown where the ionization 

35 step takes place. 
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[0031] A Fourier plane pulse shaper is preferably used with the present 
invention. Ultra-fast laser pulses contain from one to fifty optical cycles, and last only 
a few femtoseconds. This is much faster than most current electronics and therefore 
shaping with fast time gates is very difficult. On the other hand, as a consequence of 
5 the uncertainty principle, the optical spectrum spans tens to hundreds of nanometers. 
Such a large bandwidth is relatively easy to measure and to filter, and there are 
several techniques to shape the spectrum in the frequency domain, and thereby 
shape the temporal pulse upon recompression. 

[0032] In order to access the frequency domain and the individual 

10 frequency components that comprise the pulse, a geometric arrangement is 
employed, using two back-to-back spectrometers. The spectrometers are especially 
designed to introduce no net temporal dispersion: that is, all colors pass through the 
spectrometers within the same amount of time. The first spectrometer (including 
grating 25 and mirror 27) spreads the unshaped pulse spectrum along a line 

15 according to its dispersion function y(<x). The light intercepts spatial amplitude and 
phase mask pulse shaper 29 at this point. The mask output then forms the entrance 
to a second spectrometer (including grating 33 and mirror 31) which recombines the 
colors into a single shaped pulse. 

[0033] The heart of pulse shaper 29 is the programmable 256 pixel liquid- 

20 crystal mask (consisting of two overlapping 128 pixel liquid crystal arrays) that is 
placed at the Fourier plane. This mask must be capable of either attenuating the 
individual colors or shifting their phase. For alternate embodiment pulse shapers, 
two different electronically programmable masks that are capable of controlling both 
amplitude and phase have been demonstrated: a liquid crystal display ("LCD") and 

25 an acousto-optic modulator ("AOM"). A LCD pulse shaper can be obtained from CRI 
Co. and has a modulator electronic driver. 

[0034] The AOM consists of an anti-reflection coated Tellurium Dioxide 
(Te02) crystal with a piezo electric transducer glued onto one end. The central 
frequency of the acoustic wave is oc/2ti = 200 MHz. The acoustic velocity vs in the 

30 crystal is 4.2 knrvs and the light pulse spends less than 10ps in the crystal, so the 
acoustic wave moves less than 0.002 X acoustic during the transit of the light field 
through the crystal. Since the acoustic wave is essentially frozen as the optical pulse 
travels through the crystal, the complex amplitude of the acoustic wave traveling 
through the crystal in the y direction, A(t)cosoct = A(y/ys)cosccct, is mapped onto the 
35 optical field E(a) as it passes through the AOM. If some of the disoersed ontirai fteiri 
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encounters a weak acoustic wave, that frequency is attenuated; if the acoustic wave 
carrier is shifted by phase angle 0, that phase shift is imposed on the optical field. 
This pulse shaper has a total efficiency of about 20% including the diffraction 
efficiency of the AOM and the diffraction efficiency of the gratings. The diffracted 
5 light is used and the undiffracted "zero order" beam is blocked, to allow full 
modulation of both amplitude and phase in the shaped beam. The shaped beam 
than has the form 

m E shaped (Q)) = E il!pul ^)xa(co)xe i ^ ) ' 

wliere a{co)e im) = A[y(ct)) I v s \, a is the frequency, and e is a 

10 constant. 

[0035] The shaped pulses are measured using spectral interferometry. In 
this technique, the shaped laser pulse is joined with an unshaped reference pulse on 
a beam splitter, and then the combined pulses are analyzed in a spectrometer. The 
signal corresponds to a spectrally resolved interference. If the reference pulse is 

1 5 known to have a flat spectral phase, then the amplitude beating of the output beam is 
a direct measure of the spectral phase function. 

[0036] The pulse shaping devices and methods can be further achieved 
as discussed within the following publications: A.M. Weiner, "Femtosecond Pulse 
Shaping Using Spatial Light Modulators," Review of Scientific Instruments, Vol. 71, 

20 No. 5, P. 1929 (May 2000); and J.X. Tull; "High-Resolution, Ultrafast Laser Pulse 
Shaping and its Applications," Advances in Magnetic and Optical Resonance, Vol. 
20, P. 1 (1997). Fixed pulse shaping optics, such as chirped mirrors, can also be 
employed as will be discussed further hereinafter. 

[0037] . The transform limited pulses ('TL"), having all their frequencies in 

25 phase, are fed into the pulse shaper where curved mirror 27 focuses the spectrum 
onto Fourier plane 29. Changes in the phase 0 and amplitude A of the spectral 
components indicated by the computer are used to tailor the laser pulse before 
reconstruction with second curved mirror 31 and grating 33. Once compressed, the 
shaped pulse is directed to mass spectrometer 39 for evaluation. The Fourier 

30 transform relationship between the time and the frequency domain allows us to 
calculate the necessary mask to create a certain shaped pulse. These calculations 
are based on 

[2] /«-££,/«■"-* 
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and 

[3] f(t)=[ f(v)e- i2ma dv 

where v is the frequency in wave numbers, t is the time, and c is the speed of light. 
As shown in Figures 5A and B an amplitude mask blocking a portion of the spectrum 
5 of a laser pulse at the Fourier plane leads to the formation of large periodic wings in 
the shaped pulse. When the phase of different frequencies in the laser pulses 
changes in sign, as shown in Figures 6A and B, the shaped pulse becomes a pair of 
out of phase pulses. Combinations of phase and amplitude masks Can be used to 
create complex shaped pulses. The ideal pulse sequence for laser-desorption may 
10 be the combination of a long pulse to melt the matrix and vaporize the analyte 
combined with a fast pulse to cause multiphoton ionization. This sequence is based 
on the observation that below threshold laser excitation in MALDI generates a plume 
of material with very little ionization, and the fast pulse at the tail of the sequence 
would provide the prompt ionization. In Figures 7A and B, the formation of such a 
pulse is illustrated. Notice that the sharp femtosecond pulse is preceded by a long 
picosecond pulse. 

[0038] The phase and amplitude masks of the pulse shaper are controlled 
by computer. The adaptive laser source is part of a learning feedback method that 
modifies the laser pulse shape based on its success at optimizing the yield of 
20 charged proteins. Traditionally in MALDI, the laser plays a relatively passive role as 
the energy source. In the present application, the laser pulse shape takes a more 
dynamic role. Pulse shapes are envisioned which include sequences of pulses 
where each pulse in the sequence plays a different role, for example, melting, 
excitation, selective fragmentation, proton transfer and evaporation. 
25 [0039] The physical process runs itself by means of an intelligent 

"feedback" method by means of an intelligent loop. The learning method tries 
various pulse shapes, assesses their success in achieving the desired target 
excitation, and uses the knowledge gained in this way to improve the pulse shapes 
on subsequent laser shots^all with only minimal intervention of the researcher or 
system user. Changing conditions are automatically corrected within the learning 
method or feedback loop. 

[0040] Reference should now be made to Figure 4. The feedback 
software is implemented with an initial random population of pulse shapes. Each 
pulse shape is characterized by a series of numbers that specify the spectral phase 



30 
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and amplitude in each of the wavelengths or frequencies within the pulse. This 
particular parameterization for the pulse shape is itself subject to optimization; this 
makes the algorithm adaptive. The most important part of the method is the test for 
fitness of a given pulse shape. Each pulse shape is tested for its ability to generate 
5 the result that most resembles a target selected in advance. For example, the 
program will calculate the ratio between the amplitude at the desired protein signal 
and the background. Once relative success is quantified, a new generation of pulse 
shapes is produced by mating different parts of the amplitude and phase information 
from pairs of the most fit pulse shapes from the current generation. In addition, the 

10 method prescribes a small probability (5%) of random changes or mutations in the 
successful pulse shapes. Furthermore, a new set (10%) of random pulses are 
introduced in order to explore entirely new regions in the parameter space. This 
basic series of processes is iterated until the fitness converges on a "best" value. 

[0041] The convergence and robustness of the feedback method 

1 5 solutions can be measured in two different ways. FirsVthe variance in the amplitude 
and phase information itself can be monitored. As the feedback method converges 
on a solution, the values fall into a narrow range that produces the best result. Over 
the course of many generations, some of the parameters become very stable, which 
is an indication that those spectral phases and amplitudes are particularly important 

20 for driving the process that determines fitness. Secondly, the information for different 
initial conditions is monitored. If the feedback method is working it should converge 
to a similar result. 

[0042] New sets of parameters are created by modifying and combining 
elements of previous pulse shapes to create new ones. This is carried out by 

25 statistical operators that act on the phases and amplitudes of the pulse shapes. 
Operators that can be used include multi-point crossover, mutation, averaging, creep, 
smoothing, choice of phase and amplitude basis, and polynomial phase mutation. 
Crossover exchanges one or more sections of the phase and amplitude parameters 
from each of two or more pulse shapes. The resulting pulse shapes are then tested. 

30 Mutation randomly alters individual phases or amplitudes in a pulse shape. 
Averaging produces pulse "shapes by averaging the values of two or more pulse 
shapes. Creep is mutation where the final parameter value is constrained to fall 
close to the initial value. Smoothing averages nearby phase or amplitude values in 
the pulse shape. Polynomial-phase mutation produces pulse shapes by replacing a 

35 portion of the parameters with a polynomial fit. 
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[0043] A well-chosen set of operators can greatly enhance the 
performance of the feedback method and lend additional physical insight. However, 
the proper choice is usually far from obvious, so the method is allowed to adapt itself 
by letting it choose how often to use a given operator to produce new pulse shapes. 
5 The use of adaptive operators helps speed up convergence, and, perhaps more 
importantly, it helps shed light on the control mechanism at work. For example, 
crossover is more effective in the beginning of the algorithm when there is maximal 
uncertainty, since it does a good job of mixing up the information of the initial pulse 
shapes. It becomes less effective as the feedback method converges to the best 

1 0 solutions, since at this point there is much less change' in the parameters, so there is 
no longer a need to drastically change the information. Ideally, the learning program 
learns from its past mistakes and does not test possible pulse shapes which it now 
knows will fail, which saves a considerable amount of computing time. 

[0044] For each pulse shape, a number of spectra will be obtained. The 

1 5 number of laser shots that are averaged per pulse shape will depend on achievement 
of a statistically significant spectrum. At first, when the pulse shapes are the result of 
random phases and amplitudes we imagine that up to 1000 repetitions may be 
needed to distinguish the more efficient pulse shapes. This will allow one pulse 
shape per second. As the selection process proceeds large gains in efficiency can 

20 be expected. The final stages of the optimization may be carried out at a rate of 1 00 
different pulse shapes per second. The goal is to reach single pulse, femto-mol 
sensitivity. 

[0045] The learning feedback software employed in the present invention 
control system and apparatus is more fully described as follows. The preliminary 

25 investigation method and computer software steps for analyzing a pre-test unknown 
sample or specimen can be observed in Figure 9. For any new system, the test 
should start with pre-defined pulse shapes in order to obtain a basic understanding of 
the system. Among the pre-defined pulses, the shortest pulse is expected to ionize 
molecules on the surface of the sample with minimum decomposition, the longest 

30 pulse is expected to mimic the nanosecond experiments where the singly protonated 
protein may be observed. It is also noteworthy to vary the delay between two laser 
pulses from a few picoseconds to a few nanoseconds in order to appreciate the time 
scales involved. The manual inputs of steps A through C will be initially performed by 
the system operator or user through entering input data into the personal computer. 

35 The microprocessor within personal computer 37 will then control laser 23 in 
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accordance with step D, receive an essentially real time feedback input signal from 
mass spectrometer 39 in accordance with step F and then perform calculations, 
comparisons and evaluations in accordance with steps G, H and I. These automated 
steps can be substituted with manual user calculations and decisions if desired 
5 based on personal computer outputs. 

[0046] The objective of the software routine of Figure 9 is to aid in the 
selection of sample targets for further testing iterations for subsequent criteria data 
input. An optional alternate embodiment subroutine includes shooting long laser 
beam pulses then quick short laser beam pulses, with a separation set by an optical 

10 delay of less than ten nanoseconds. The short pulse of approximately 50 
femtoseconds is performed in order to look for fragmehtation and the matrix mass. 
Laser beam pulses of between approximately ten and 100 picoseconds are 
performed to look for the parent mass. The ultrafast laser beam pulse durations 
employed with the present invention advantageously allow for approximately 1000 

15 laser beam shots at a single sample or specimen without significant degradation of 
the specimen; this allows for quicker and less expensive usage of the apparatus 
while also encouraging statistically more accurate results. The long and short pulse 
combinations can be used in addition to or without the benefit of pulse shaping. 
. Otherwise, the control system and apparatus are the same as discussed herein. 

20 [0047] The real time learning feedback method and computer software 

will now be described in greater detail with regard to Figures 10A and 10B. This 
method and software are employed to statistically optimize the repetitive identification 
of molecularly complex and unknown samples or specimens in a highly automated 
and relatively quick manner. The data that is obtained in the mass spectrometer for 

25 each laser pulse consists of a two-dimensional array of numbers which are stored in 
. random access memory of the personal computer. The first column of the array 
contains arrival times for the ions, when the data is obtained from the time-of-flight 
mass spectrometer. Equivalent numbers can be obtained from different mass 
spectrometry units such as quadrupole or ion-cyclotron spectrometers. This 

30 numbers can be converted to a mass assuming a charge for the species. The 
second number corresponds to the intensity (Rl) of signal at that particular arrival 
time. This number correlates with the probability of detecting a charge species with 
that mass (see, for example, Figure 8). The data acquisition involves collecting a 
number of data sets and then calculating the sum; this action is performed by the 

35 personal computer. The target is an ideal mass spectrum defined by the system 
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user. The target can be the maximization or minimization of a signal at a particular 
single mass, or a number of masses with or without background suppression. 

[0048] Each pulse shape is defined by a two-dimensional array of 
numbers which are stored in random access memory of the personal computer or 
5 otherwise accessible on a read only basis by the personal computer from a disk or 
the like. The length of each column determines the resolution with which the 
spectrum of the laser pulse is being sculpted. For a liquid spatial-light modulator, the 
number of pixels typically determines this number. The first column of numbers 
determines the amplitude transmission coefficient for every pixel. The second 

10 column determines the phase delay for every pixel. The entire array of amplitudes 
and phases determines the final shape of the output pulse. The shortest pulse is 
pre-defined as the shortest duration possible for the laser system will the longest 
pulse is pre-defined as the longest pulse that can be made with the pulse shaper. A 
two-pulse combination is pre-defined as the combination of an unshaped pulse with a 

15 shaped pulse. Pre-defined ultraviolet or infrared pulses at 400 or 800 nm, for 
example, can be used. 

[0049] For the first iteration, the computer generates a number of different 
seed pulsed shapes with a pre-determined resolution using a random number 
generator. For each pulse shape, a number of mass spectra (the number of 

20 averages) are summed. The number of averages is determined by the reproducibility 
of the data obtained. The resulting average data is compared with the target and a 
success value is calculated] It is expected as the optimum pulse shape is achieved, 
the number of averages can be reduced because of the higher efficiency of these 
pulses. A number of successor pulses are chosen based on their success value and 

25 the rest of the pulse shapes are discarded. The best one is stored as the optimum. 
The number of "more iterations" is reset to zero; this number keeps track of how 
often the optimum pulse is redefined and is used to determine convergence. The 
success value is a number between zero and one that is assigned to each average 
set; this number quantifies how well the average data set approximates the target 

30 and can be obtained, for example, from a sum of the differences or a sum of the 
differences squared. Step A, inputting/defining the target value based on the mass, 
the percent of mass and minimum ion count level is very important as one wants to 
define regions in the spectrum to be maximized and/or minimized. This also assists 
in reducing background noise or interference from the mass spectrometer detections. 

35 Step B, inputting/defining the number of shots to average for each seat pulse can be 
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set, for example, at 1000 laser beam shots for a single sample without significant 
degradation; this is highly advantageous over the typical 100 or less maximum shots 
traditionally used with prior nanosecond lasers before the specimen is unusable. 

[0050] The adjustment and focusing of the shaped beam onto the sample 
5 in step J, can be performed manually or automatically by the personal computer. For 
example, a relatively weak diode laser, having the same wavelength as the main 
femtosecond laser and following the same path, can be used with a CCD camera to 
aim and focus the main shaped laser beam onto the sample in an automated and 
computer controlled, real time feedback manner. Furthermore, step L allows for a 
1 0 setting of delayed pulse extraction within the mass spectrometer. The calculation, 
comparison, and determination steps, such as those of steps N through V, are all 
conducted in an automatic manner within the microprocessor of the personal 
computer. 

[0051] Step R allows the microprocessor to determine the best pulse 

15 shape with the highest success value and store it as the optimum value in the 
random access memory of the central processing unit. The computer will then pick 
approximately the ten percent best pulse shapes based on the highest success 
values and then reset the discarded values and automatically generate new laser 
pulse shapes in step V for subsequent testing iterations on the same specimen. The 

20 generation of new seed pulses is important for the success of the feedback method. 
The goal is to arrive at the optimum pulse in the shortest number of iterations making 
sure that one has searched the entire range of parameters, the global maximum. 
The "cost functional" refers to the statistical pressure that is placed on the optimum 
pulse shape in order to simplify it. For example, once an optical pulse shape or other 

25 characteristic is found, it may be important to determine how sensitive the outcome is 
to each of its amplitude and phase components. Perhaps a much lower resolution 
can produce the results. The simpler the pulse the easier it is to reproduce and 
interpret the results in terms of physical concepts. For certain cases, the shape can 
be simple enough that it can be prepared without a. shaper which would allow for a 

30 less expensive alternative to the preferred pulse shaping; for example, a combination 
of two or three different 800 nm pulses, or a combination of infra-red and ultra-violet 
pulses could be employed as a modified or optimized pulse. Once statistical 
convergence has been determined by the personal computer, then the test is 
complete by determining the optimum pulse characteristics (whether they be pulse 
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shape, pulse duration or any other such variable laser beam characteristic) for the 
corresponding and now post-test identified specimen. 

[0052] The time scale of some of the processes that occur during MALDI 
may be longer than the femtosecond pulses. In a first variation, the pulse shaper can 
5 be used to produce pulse sequences up to ten picoseconds apart. Optical delay 
lines can be used to increase this time delay in the nanosecond range if needed. In 
a second variation, the wavelength of the pulses being shaped is 800 nm. A second 
harmonic crystal is all that is needed to convert the wavelength to 400 nm, however, 
the shaper is capable of regulating the energy delivered to the sample without 
1 0 changing the carrier frequency (wavelength) of the laser. 

[0053] The sensitivity and flexibility gained should make this unit cost 
effective, especially if in-source selective bond cleavage is achieved. Another 
variation, used to further minimize cost, provides that optimal pulse shapes may be 
synthesized from a combination of less expensive laser sources. 
15 [0054] The concern of missing the MALDI crystal by using a' sharply 

focused laser, usually observed when a weak and inexpensive laser source is used, 
can be overcome with the present system since it is able to use focal spots as large 
as a millimeter. The peak intensity of the laser will exceed 10 11 W/cm 2 , the ionization 
threshold due to multiphoton excitation, using 0.1 mJ per pulse focused on a 1 mm 
20 diameter spot. 

Protein Sequencing 

[0055] Laser desorption mass spectrometry can be employed with the 
present invention for identification and protein sequencing. This is significantly 
enhanced and made possible by the ultra-fast laser pulses and learning feedback 
25 system used. The matrix has been shown to enhance the yield of charged protein for 
analysis by MS detection. The matrix: phosphor diester backbone interaction has 
been shown to play an important role. The use of liquid matrices such as glycerol 
and lactic acid for IR-MALDI may bring some additional flexibility to sample 
preparation and delivery to the MALDI instrument. The "Ladder Sequencing" method 
30 involves a partial Edman degradation with phenyl isothiocyannate and using phenyl 
isoccyanate as a terminating agent. Partial enzymatic hydrolysis of polypeptides 
using trypsin is another strategy for protein sequencing. Trypsin digestion attaches 
only bonds in which the carboxyl group is contributed by either a lysine or an arginine 
residue. Analysis of metastable species in MALDI-PSD using a reflection TOF 
spectrometer leads to valuable structural information. The introduction of 'delayed 
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extraction' in MALDI allows improved resolution, suppression of matrix background, 
reduction of chemical noise, and minimization of the effect of laser intensity on 
performance. 

[0056] MALDI is a soft ionization technique which produces protonated 
5 molecules that undergo very little or no subsequent fragmentation due to the low 
amount of energy imparted during the ionization process. Therefore, MALDI can be 
used to analyze mixtures of peptides because the mass spectrum of one peptide is 
unlikely to overlap with the spectrum of another. Ideally, cleavage of the ionized 
peptide at each peptide bond would provide a mass spectrum that could be 

10 interpreted, using knowledge of the masses of the amino acid residues, to deduce 
the sequence. However, as conceptually illustrated in Figure 11, cleavage on either 
side of the I-carbon is also possible to give fragment ions, which, while diagnostically 
useful, also complicate the spectrum. It is also noteworthy that cleavage at any 
designated bond can generate either an N-terminal ion (a, b, c) or a C-terminal ion (x, 

15 y, z), the predominance of which for a protonated peptide (MH + ) depends on the 
locus of the more basic residues. In reality, the fragmentation process is more 
complicated than suggested in Figure 1 1 ; for example, creation of a y-ion involves 
hydrogen transfer from the N-terminal side of the peptide bond and retention of the 
ionizing proton. In addition, there can be fragmentation of the side chain on certain 

20 residues; for example, fragmentation involving cleavage at the ^carbon of leucine 
and isoleucine generates w-ions, which distinguish these two isomeric residues. 

[0057] Recognizing the ion types as represented by the appearance of 
peaks in the mass spectrum is not critical, as most strategies for interpretation, 
especially those using an algorithm, involve an iterative computational approach. 

25 However, the beginning of a C-terminal series of fragments can be distinguished 
from the start of an N-terminal series. The largest b-ion will be represented by a 
peak at high m/z value that differs from that representing MH + by a number of mass 
units equal to the sum of the mass of an amino acid residue plus the mass of water 
due to expulsion of the C-terminal residue, which contains the hydroxyl group. On 

30 the other hand, the largest y-ion is represented at a high m/z value by a peak 
differing from that for MH + by a number of mass units equal to only the mass of an 
amino acid residue. 

[0058] In principle, the sequence of a peptide is deduced from a mass 
spectrum in which a complete series of any given ion type are represented. In 
35 practice, however, a complete series of any one type is rarely observed, but in 
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96. The communications system of Claim 93 wherein the laser is a 
femtosecond laser and the first device includes a chirped phase mask element with a 
passive pulse changing characteristic and a dispersive element. 

97. A method of controlling a laser beam, the method comprising: 
(a) emitting a laser beam pulse of less than about 1 picosecond 



(b) shaping the laser beam pulse with a smooth phase function; 

(c) controlling the nonlinear optical processes induced by the laser 



duration; 
and 

10 beam pulse. 

98. A method of ionizing and determining a characteristic of a 
specimen, the method comprising: 

(a) emitting a laser beam having a pulse duration of less than 
about 51 femtoseconds; 

15 (b) shaping the laser beam pulse; 

(c) sensing a characteristic of the specimen after the laser beam 
pulse has ionized at least a portion of the specimen; 

(d) varying laser beam pulse shapes emitted upon the specimen 
for subsequent iterations; 

20 ( e ) . sensing characteristics of the specimen after the laser beam 

has ionized at least a portion of the specimen for the subsequent iterations; and 
(f) comparing the sensed results. 

99. The method of Claim 98 further comprising determining a 
protein sequence of the specimen based at least in part on the sensed 

25 characteristics.. 
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